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A summary is provided for the First AIAA Sonic Boom Workshop held 11 January 
2014 in conjunction with AIAA SciTech 2014. Near-field pressure signatures extracted 
from computational fluid dynamics solutions are gathered from nineteen participants rep- 
resenting three countries for the two required cases, an axisymmetric body and simple delta 
wing body. Structured multiblock, unstructured mixed-element, unstructured tetrahedral, 
overset, and Cartesian cut-cell methods are used by the participants. Participants provided 
signatures computed on participant generated and solution adapted grids. Signatures are 
also provided for a series of uniformly refined workshop provided grids. These submissions 
are propagated to the ground and loudness measures are computed. This allows the grid 
convergence of a loudness measure and a validation metric (difference norm between com- 
puted and wind tunnel measured near-field signatures) to be studied for the first time. 
Statistical analysis is also presented for these measures. An optional configuration includes 
fuselage, wing, tail, flow-through nacelles, and blade sting. This full configuration exhibits 
more variation in eleven submissions than the sixty submissions provided for each required 
case. Recommendations are provided for potential improvements to the analysis methods 
and a possible subsequent workshop. 
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Wing span 

Freestream pressure difference normalized by freestream 

Characteristic spacing of the grid 

Distance below the model 

Coverage factor 

Model reference length 

Pressure 

Radius from the model centerline 
Validation metric 

Forward validation metric integration extent 
Aft validation metric integration extent 
Tangent of freestream shock angle 
Estimate of sample median 
Estimate of sample standard deviation 
X , shifted and normalized by ( X — /3R)/L 
Off-track angle (zero is directly below model) 

Freestream conditions 
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I. Introduction 


The First American Institute of Aeronautics and Astronautics (AIAA) Sonic Boom Prediction Workshop 
(SBPW-1) was held on January 11, 2014 at National Harbor, Maryland. There were approximately 50 
attendees from 7 nations. The objective of this workshop was to assess the state of the art for predicting 
near-field pressure signatures needed for accurate and reliable sonic boom prediction. These near-field 
pressure signatures, extracted from Computational Fluid Dynamics (CFD) solutions, were gathered from 
international participants. The participants utilized a series of uniformly refined workshop provided grids 
and grids created with each participant’s best practices for computing solutions on the provided geometries. 

Statistical analysis of the SBPW-1 submissions is facilitated by reducing the submitted near- field sig- 
natures to a set of scalars for N-version testing. 1 SBPW-1 focuses on near- field CFD, but predicting the 
acoustic signature on the ground and how it is perceived by humans is the goal of the sonic boom community. 
The submitted near- field signatures are propagated to the ground (Fig. 1) and noise measures are computed 
to reduce these signatures to relevant scalars. The comparison of CFD to wind tunnel measurements is 
also of interest to the sonic boom community, so a validation metric is computed to quantify the difference 
between near- field CFD and measurements. The workshop requirement of computing solutions on uniformly 
refined grids allows, for the first time, the grid convergence of these scalars to be examined. The statistical 
products and grid convergence studies begin the process of quantifying the uncertainty of near-field CFD 
employed by the international participants. 


CFD Domain 



Waithe 2 3 introduced the workshop in special session at the 51 st AIAA Aerospace Sciences Meeting, a 
year before the workshop. The name of the workshop was changed and one of the required cases was 
replaced with the SEEB-ALR (computer program SEEB 4 modified with Lockheed-Martin proprietary Aft 
Lift Relaxation) because as-built geometry was available and the model was larger in size. The larger model 
had longer shaped signature that enabled more accurate wind tunnel measurements. The DWB (69° Delta 
Wing Body) 5 was also chosen as a required case and the LM1021 (Lockheed-Martin 1021-01 configuration) 6 
was an optional case. See Fig. 2 for an illustration of the models. This special session included a summary 
of the 2008 NASA Sonic Boom Prediction Workshop, 7 which also utilized the DWB model at a different 
flow condition, a number of axisymmetric bodies, and a configuration with wing, body, tails, and nacelles. 
The 2008 workshop had a much smaller scope than SBPW-1; it was limited to five NASA participants and 
the only analysis of results was comparing near-body signatures to wind tunnel measurements. SBPW-1 
benefits from the lessons learned at the 2008 NASA Sonic Boom Prediction Workshop and previous AIAA 
workshops. 8-11 

The data and analysis methods of this summary are intended to contribute toward the discussion of 
replacing the prohibition of overland supersonic flight 13 with a certification standard. The existing ban on 
overland flight increases the cost and block time 14 of supersonic fights between city pairs and is an important 
driver of the economic viability of a supersonic transport. 15,16 
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II. Evaluation Methods 

The SBPW-1 participants submitted near-field pressure signatures: normalized delta pressure as a func- 
tion of x. While these signatures could be compared directly to each other and wind tunnel measurements, it 
was helpful to reduce the properties of these signatures to scalars for grid refinement studies and application 
of statistical methods. This data reduction was performed by the authors, not by the participants, in an 
attempt to apply these methods as uniformly as possible to the submissions. 

II. A. Loudness 

Loudness and annoyance are inherently subjective measures of a sonic boom experience. Many noise de- 
scriptors have been evaluated as loudness predictors in human experiments, for example see Leatherwood 
et al. 17 This summary utilizes two of them. The Mark VII perceived loudness level (PL) of Stevens 18 is 
adopted as a measure to compare submissions because it is correlated with loudness and annoyance in mul- 
tiple experiments. 17 The C-weighted sound exposure level (CSEL) is also used because of its emphasis on a 
lower and wider range of frequencies. The method of Shepherd and Sullivan 19 is used to calculate PL and 
CSEL from ground signatures. The propagation of the near-field submissions to the ground is described in 
the next section. 

II. B. Propagation 

Assumptions are made to scale the wind tunnel models and test conditions to a notional full-scale vehicle for 
propagation and loudness calculations. The wind tunnel test Mach number is used for the full-scale vehicle. 
All full-scale vehicles are assumed to be at 55,000 foot altitude in a no- wind US Standard Atmosphere. 20 The 
1.6 Mach SEEB-ALR wind tunnel model is 0.006 scale, 21 the 1.7 Mach DWB wind tunnel model is 0.0065 
scale, 3 and the 1.6 Mach LM1021 is 0.008 scale. 21 Before propagation, the near-field pressure signature is 
scaled and the initial and final delta pressures are linearly ramped to zero, windowing the submission and 
wind tunnel measurements. The CFD submission computed with the largest grid for each configuration is 
shown with the wind tunnel (WT) measurement at the same location in Fig. 3. This windowing forces the 
initial and final pressures to freestream in an attempt to reduce the impact of wind tunnel mounting hardware 
on propagated signatures. Participants extracted pressure between different forward and aft locations. This 
windowing also eliminates these differences between submissions. 

The sBOOM code of Rallabhandi 22 is used to propagate near-field signatures to the ground using an 
augmented Burgers equation solver. This is the atmospheric propagation step in Fig. 1. Rallabhandi 23 
provides additional details of sBOOM including an adjoint formulation. 23 Version 1.4 of sBOOM is used 
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Figure 3. Full-scale near-field signature with windowing. 
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in this summary. It utilizes the atmospheric humidity model from ANSI SI. 26, Annex C. 24 The previous 
version 1.3 of sBOOM was used for the workshop overview and summary presentations. It had a bug in the 
atmospheric pressure calculation based on the hydrostatic equation and a different atmospheric humidity 
model. The changes for sBOOM version 1.4 resulted in a median PL increase that is approximately 1.2 dB 
louder for the SEEB-ALR and approximately 0.8 dB louder for the DWB than the statistics computed with 
version 1.3. 

Like the results of any finite- difference method, the ground signatures predicted by sBOOM has a sen- 
sitivity to the number of points used to discretize the signature. 22 For example, Cleveland 25 shows the 
dependency of the finite- difference molecular relaxation terms on the grid resolution. To verify the consis- 
tency of the sBOOM finite difference method, a series of uniformly refined grids are used. This refinement 
can be expressed as spatial refinement in the aircraft’s reference frame (points per body length) or temporal 
refinement (sampling frequency) for a stationary (ground) observer. The resolution of the grid is expressed 
in points per body length and frequency in Fig. 4. Spatial resolution is important to the near-field CFD 
pressure extraction process and temporal resolution is important to interpreting results in the frequency 
space of the loudness calculation. PL as a function of sBOOM grid refinement is shown for wind tunnel 
measurements (WT) and the participant signature computed on the largest grid (CFD). The PL based on 
the LM1021 wind tunnel configuration is louder than predictions of a full-scale free-flight model due to a 
number of factors; see section V.C for details. 

The sBOOM grid has 36,000 points for the generation of participant statistics, which corresponds to 
18,000 points per body length (90K Hz sampling frequency) for the SEEB-ALR, 12,000 points per body 
length (161K Hz sampling frequency) for DWB, and 18,000 points per body length (100K Hz sampling 
frequency) for the LM1021. At this resolution, the computed PL is within a tenth of a dB of the asymptotic 
PL value. The LM1021 shows slightly different convergence behavior for CFD and WT below 18,000 points 
per body length (100K Hz sampling frequency). There is a nearly constant difference between the CFD and 
WT signatures for a range of resolutions for the SEEB-ALR and DWB, which indicates the statistics might 
be insensitive to resolution within this frequency range. 

The default action of sBOOM is to detect the initial shock and interpolate the near field signature aft of 
this first shock to the available points. This was disabled to ensure that all signatures have equal resolution, 
because undershoots or smeared front shocks would result in a different effective resolution (points per body 
length). 

The frequency spectra of the CFD ground signatures is shown in Fig. 5 when propagated with two 
different resolutions. This resolution is 18,000 points per body length (31 KHz) and 6,000 points per body 
length (93 KHz) for the SEEB-ALR and 12,000 points per body length (161 KHz) and 4,000 points per body 
length (53 KHz) for the DWB. The finer sBOOM grid only affects the higher frequencies (above 1 KHz). 
The SEEB-ALR has more energy in the lower frequencies (below 10 Hz) than the DWB. The DWB has more 
energy above 10 Hz than the SEEB-ALR. The differences in the spectra are due to the SEEB-ALR flattop 
signature and the multiple N-wave character of the DWB (shown in Fig. 3), which agree with the analytic 
predictions of Shepherd and Sullivan. 19 

II. C. Validation Metric 

A validation metric is employed to quantify the difference between an experimental measurement and a 
computational prediction. This technique was also employed in the analysis of the AIAA Shock Boundary- 
Layer Interaction Workshop submissions. 11 For a detailed discussion of the construction and use of validation 
metrics, see Oberkampf and Roy. 26 The wind tunnel measurements in this summary are the mean of an 
ensemble of measurements. See section VI for details of this averaging process used to extract the near- 
field signature from the noisy environment in supersonic wind tunnel test section. The mean wind tunnel 
measurements are shifted in x so that the location of the SEEB-ALR and LM1021 nose shocks and the 
DWB main expansion of the measurement mean matches the median of the participants. This difference is 
evaluated as an integral, 


V = 



Ap\ 

Poe) 


submission 


Ap\ 

Poo J measurement 


(1) 


The xq and x\ locations define the same window used for each configuration during scaling for atmospheric 
propagation. Numerically, this integral is evaluated with the trapezoidal rule after all of the data points from 
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Figure 4. Perceived Level as a function of sBOOM grid resolution. 




(a) SEEB-ALR. (b) DWB. 

Figure 5. Ground signature frequency spectra at two sBOOM grid resolutions. 
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the submission and mean wind tunnel measurement are linearly interpolated to the common x locations of 
the submission and mean wind tunnel data. 

III. Statistical Approach 

The intention of this applied statistical approach is to measure the repeatability of an experiment that 
is run multiple times, where a workshop submission is treated as an experiment. In this context, no sub- 
mission is treated as the “correct” prediction because all submissions contains known and unknown sources 
of uncertainty. The approach is based on the principles of the AIAA Drag Prediction Workshop (DPW) 
statistical summary. 1,9 Two samples of the submissions are examined: a sample where the assumption of 
independence is violated by including uniformly refined grids from the same code and participant to increase 
the sample size, and a sample where the coarser uniformly refined grids are excluded. For each population, 
the scatter limits are established as 

(i d= Ka. (2) 

The median is used to estimate the sample mean (i because it is robust to outliers. The sample standard 
deviation a is estimated with the corrected sample standard deviation. A coverage factor of K — \f?> is 
calculated by assuming a uniform distribution of results because the actual distribution of the submissions 
is unknown. 9 The first two DPW statistical summaries assumed 99% coverage of a normal distribution, 
which resulted in a larger K. 1 A submission is significant and considered outlier if it is outside the limits 
established by Eq. (2). A larger choice of K would result in fewer submissions being outliers and a smaller 
choice of K would result in more submissions being outliers. 

IV. Participants 

The groups contributing to the workshop are listed in Table 1. The group letter identifier is included 
in the submission identifier throughout this summary. Some groups contain a single contributor, but most 
contain two or more collaborators. The groups have members from the United States of America, France, 
and Japan. Government agencies, industry, and academia are represented. Each group submitted results and 
prepared a presentation for the workshop. All group presentations are available on an FTP server, ftp:// 
lbpw-ftp.larc.nasa.gov/outgoing/lbpwl/presentations/, with the workshop introduction, overview, 
and summary presentations. 


Table 1. Participant groups. 


Group 

Organizations 

Contributors 

A 

NASA Langley 

M. Park 

B 

NASA Langley 

S. Rallabhandi 

C* 

NASA Langley, NASA Ames 

A. Elmiligui, S. Cliff 

D 

INRIA 

A. Loseille, F. Alauzet 

E 

Dassault 

F. Dagrau 

F 

ONERA 

1. Salah El Din, A. Minelli, R. Grenon, G. Carrier 

G 

University of Miami 

J. Gan, G. Zha 

H 

Stanford University, ONERA, INRIA 

F. Palacios, T. Lukaczyk, J. Alonso, G. Carrier, A. Loseille 

1 

Lockheed Martin 

J. Morgenstern 

J 

JAXA 

H. Ishikawa, Y. Makino 

K 

Tokyo University of Science, JAXA 

Z. Lei, Y. Makino 

L 

JAXA 

M. Kanamori, A. Hashimoto, T. Aoyama 

M 

NASA Ames 

M. Aftosmis, M. Nemec 

N 

Boeing 

T. Magee, S. Fugal, D. Lazzara, M. Sexton, E. Unger 

0 

Gulfstream 

D. Howe 

P 

NASA Langley 

W. Li 

Q 

Lockheed Martin 

M. Buonanno, F. Marconi 

R 

NASA Ames 

J. Housman, E. Sozer, C. Kiris 

S* 

Arizona State University 

J. Jensen 


Group S submitted results after the workshop, but was included in the Group C workshop 
presentation. 
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V. Models and Required Grids 


The models for SBPW-1 were chosen to provide a range of geometric complexities, see Fig. 2. Grids 
were provided for both the required and optional cases in order to gain the largest degree of participation. 
Two cases were required: the SEEB-ALR and DWB. Participants were required to run at least one of the 
provided grids for these cases. The provided unstructured mixed-element and tetrahedral-only grids for the 
SEEB-ALR and DWB were described by Park et al. 27 A series of uniformly refined grids were provided to 
quantify the spatial convergence of the near-field CFD schemes. The wide range of grid resolutions were 
practical for both node-centered and cell-centered schemes. Care was taken to ensure that the relative grid 
spacing distribution remained as similar as possible for the unstructured family of grids. The number of 
nodes and elements approximately doubled and a characteristic grid spacing decreased by the cube root of 
one half between successive grids in the series. 27 

Structured grids were also provided with instructions to recursively remove every other point in each 
computational direction to produce three coarser grids. Participants were encouraged to provide their own 
grids generated with best practices or solution adaptive methods. A single tetrahedral LM1021 grid was 
provided for the optional case as described by Morgenstern, Buonanno, and Marconi. 6 

V.A. SEEB-ALR 

The SEEB-ALR model is described by Morgenstern et al. 21 It is designed as an axisymmetric body, but 
slight asymmetry and other imperfections were introduced during wind tunnel model construction. These 
small imperfections in the geometry can be seen as high frequency, small amplitude pressure fluctuations 
in the flattop portion of the near-field signature, see Fig. 3(a). Resolving these pressure fluctuations is a 
deceptively difficult task and makes this an interesting configuration for the workshop. The SEEB-ALR solid 
model is constructed in NX® from a STEP file of the scanned, as-built wind tunnel model. The geometry 
differences between the as-designed and as-built are summarized by Cliff et al. 28 and Aftosmis and Nemec. 12 
The model has a reference length of L = 17.68 in. 

The participants provided submissions at H — 21.2 in. and H = 42.0 in. below the centerline of the 
model, see Table 2. The SEEB submissions are listed in Table 3. The Groups are listed in Table 1. The 
Method is the CFD solver employed and Physics denotes the linearized Euler, Euler, Laminar, or Reynolds 
Averaged Navier-Stokes (RANS) with the specified turbulence model. Euler results were requested by the 
committee and the committee provided grids had a normal spacing intended for Euler calculation, so that 
is the most common submission. Type is a description of the construction of the grid. The Grid identifier 
is the grid generation (or adaptation) method: Output is output-based (goal-based) adaptation, PW is 
Pointwise™ commercial grid generation tool, SI is the finest provided multiblock structured grids, and S2, 
S3, and S4 are recursively refined. The grids starting with M are committee-provided 27 mixed-element grids 
and the grids starting with T are these mixed-element grids converted to all tetrahedra. The three numbers 
after the letter indicate the characteristic grid spacing h, where a smaller h has more control volumes. The 
final column is the number of control volumes, which may be the number of nodes or cells depending on the 
CFD scheme. The submissions are provided for a large range for grid sizes. The smallest has 14 thousand 
control volumes and the largest has 73 million. 

Table 2. SEEB-ALR pressure signature extraction locations. 

H H/L 0 Available Wind Tunnel Measurements 

21.2 in. 1.2 0° Yes 

42.0 in. 2.4 0° Yes (poor quality) 


V.B. DWB 

The DWB (Fig. 6) was originally described as Model 4 in Hunton, Hicks, and Mendoza. 29 It was selected 
for the workshop because it has a simple axisymmetric fuselage and diamond airfoil wing with available 
wind tunnel data both on and off the centerline. It provides a slight variation in on- and off-track near-field 
signature because the wing has different impact on the signature at different 0 angles. This model was also 
of interest because many researchers have used it as a validation case, and the 2008 NASA Sonic Boom 
Prediction Workshop examined this model at angle of attack five years before SBPW-1. 7 
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Table 3. SEEB-ALR submissions. 


Submission 

Group 

Method 

Physics 

Type 

Grid 

Cont. Vol. 

1 

A 

CFL3D 

Euler 

Struct 

SI 

7,323,648 

2 

A 

CFL3D 

Euler 

Struct 

S2 

915,456 

3 

A 

CFL3D 

Euler 

Struct 

S3 

114,432 

4 

A 

CFL3D 

Euler 

Struct 

S4 

14,304 

5 

A 

FUN3D 

Euler 

Mixed 

M080 

10,159,421 

6 

A 

FUN3D 

Euler 

Mixed 

M100 

5,252,466 

7 

A 

FUN3D 

Euler 

Mixed 

M125 

2,537,977 

8 

A 

FUN3D 

Euler 

Mixed 

M156 

1,165,851 

9 

A 

FUN3D 

Euler 

Mixed 

M200 

674,734 

10 

A 

FUN3D 

Euler 

Tet 

T080 

10,159,421 

11 

A 

FUN3D 

Euler 

Tet 

T100 

5,252,466 

12 

A 

FUN3D 

Euler 

Tet 

T125 

2,537,977 

13 

A 

FUN3D 

Euler 

Tet 

T156 

1,165,851 

14 

A 

FUN3D 

Euler 

Tet 

T200 

674,734 

15 

B 

FUN3D 

Euler 

Tet 

VGRID 

21,519,091 

16* 

C 

USM3D 

Euler 

Tet 

T080 

59,565,426 

17* 

C 

USM3D 

Euler 

Tet 

T100 

30,647,394 

18* 

C 

USM3D 

Euler 

Tet 

T125 

14,695,383 

19* 

C 

USM3D 

Euler 

Tet 

T156 

6,662,284 

20* 

C 

USM3D 

Euler 

Tet 

T200 

3,840,249 

21 

C 

USM3D 

Euler 

Tet 

PW/MCAP 

2,236,059 

22 

C 

USM3D 

Laminar 

Tet 

PW/MCAP 

6,196,118 

23 

C 

USM3D 

SA 

Tet 

PW/MCAP 

6,196,118 

24* 

D 

WOLF 

Euler 

Tet 

Output 

208,280 

25* 

D 

WOLF 

Euler 

Tet 

Adapted 

875,511 

26* 

D 

WOLF 

Euler 

Tet 

T080 

10,159,421 

27* 

D 

WOLF 

Euler 

Tet 

T100 

5,252,466 

28* 

D 

WOLF 

Euler 

Tet 

T156 

1,165,851 

29* 

D 

WOLF 

Euler 

Tet 

T200 

6,747,34 

30 

E 

Eugenie 

Euler 

Tet 

T080 

10,159,421 

31 

E 

Eugenie 

Euler 

Tet 

T100 

5,252,466 

32 

E 

Eugenie 

Euler 

Tet 

T125 

2,537,977 

33 

E 

Eugenie 

Euler 

Tet 

T156 

1,165,851 

34 

E 

Eugenie 

Euler 

Tet 

T200 

674,734 

35 

F 

CEDRE 

Euler 

Tet 

T100 

30,647,394 

36 

F 

CEDRE 

Euler 

Tet 

T156 

6,662,284 

37 

F 

CEDRE 

Euler 

Tet 

T200 

3,840,249 

38 

G 

FAS IP 

Euler 

Struct 

SI 

7,323,648 

39 

H 

SU2 

Euler 

Mixed 

M080 

10,159,421 

40 

H 

SU2 

Euler 

Mixed 

M100 

5,252,466 

41 

H 

SU2 

Euler 

Mixed 

M156 

1,165,851 

42 

H 

SU2 

Euler 

Mixed 

M200 

674,734 

43 

1 

Linear 

Linear 

N/A 

N/A 

N/A 

44 

J 

TAS 

Euler 

Mixed 

M080 

27,865,176 

45 

J 

TAS 

Euler 

Mixed 

M100 

14,325,194 

46 

J 

TAS 

Euler 

Mixed 

M200 

1,793,549 

47 

J 

UPACS 

Euler 

Struct 

SI 

7,323,648 

48 

J 

UPACS 

Euler 

Overset 

Fixed 

27,646,000 

49 

K 

ADCS 

Euler 

Struct 

SI 

7,986,107 

50 

L 

FaSTAR 

Euler 

Mixed 

M100 

10,159,421 

51 

L 

FaSTAR 

Euler 

Mixed 

M125 

5,252,466 

52 

L 

FaSTAR 

Euler 

Mixed 

M156 

1,165,851 

53 

L 

FaSTAR 

Euler 

Mixed 

M200 

674,734 

54 

M 

Cart3D 

Euler 

Cart 

Output 

2,012,184 

55 

N 

Cart3D 

Euler 

Cart 

Fixed 

11,147,363 

56 

N 

OVERFLOW 

Euler 

Overset 

PW 

26,309,760 

57 

O 

CartOver 

Euler 

Hybrid 

Fixed 

73,000,000 

58 

P 

Cart3D 

Euler 

Cart 

Fixed 

18,096,942 

59 

Q 

CFD++ 

SA 

Struct 

Gridgen 

2,000,000 

60 

R 

LAVA 

Euler 

Overset 

Fixed 

21,800,000 

61 

R 

LAVA 

SA 

Overset 

Fixed 

21,800,000 

62 

R 

LAVA 

Euler 

Struct 

Axi 

330,000 

63 

R 

LAVA 

SST 

Struct 

Axi 

330,000 

64* 

S 

OVERFLOW 

SA 

Overset 

Fixed 

21,300,000 


* Updated or submitted after the workshop. 
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(a) DWB planform, dimensions in cm. 29 


(b) DWB wind tunnel model and sting. 


Figure 6. Delta Wing Body and wind tunnel mounting sting. 


The DWB fuselage has an aft-facing step base region where it attaches to the wind tunnel sting, see 
Fig. 6(b). The solid model of the DWB is extracted from a STEP file exported from CATIA® V5, where 
it was created from the analytic description in the original report. 5 The sting geometry was created from 
drawings used to machine the wind tunnel sting. A simplification was made by omitting the strain gauges, 
wires, and instrumentation slots visible in a photograph of the model, see Fig. 6(b). The sting was inserted 
0.079184 meters (3.1175 inches) into the wing body and a Boolean addition produced the solid model. The 
model has a reference length of L = 0.1752 m and a half span of 6/2 = 0.0345 m, which produces a reference 
length to half span ratio of 5.1. 

The participants provided submissions containing ten extraction locations both on and off the centerline 
of the model, Table 4. The DWB submissions are listed in Table 5. See section V.A for a description of the 
abbreviations that are common for the SEEB-ALR and DWB tables. 

Table 4. DWB pressure signature extraction locations. 


H 

H/L 

H/(b/ 2) 

<t> 

Available Wind Tunnel Measurements 

0.5 in. 

0.07 

0.34 

0° 

No 

21.2 in. 

3.07 

13.4 

0° 

No 

24.8 in. 

3.60 

16.7 

0°, 30°, 60°, 90° 

Yes 

31.8 in. 

4.61 

21.5 

0°, 30°, 60°, 90° 

Yes 


V.C. LM1021 

The NASA High Speed Project sponsored wind tunnel testing to validate very low sonic boom designs. 21 
The resulting measurements were found to match CFD predictions within 1 PL (dB) when propagated to the 
ground. 6 Based on this success, the LM1021 model (Fig. 7) was selected as a case for the workshop. It was 
a challenging candidate because analysis performed during model development indicated that the near-field 
signature changed when laminar or turbulent CFD was applied. This sensitivity to boundary layer state 
implied that near-body signatures extracted from Euler calculations could be different from viscous CFD 
calculations. Extensive wind tunnel model oil surface flow photos were captured for comparison with CFD 
analyses and to provide an understanding of boundary layer state and flow features. This full configuration 
model included wing, body, V-tail, and three nacelles that provide a signature that varies with off-track 
angle. 




Figure 7. Perspective and three-view of the LM1021 free-flight model. 

For non-axisymmetric bodies, the spanwise distribution of acoustic disturbance sources reach a near-field 
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Table 5. DWB submissions. 


Submission 

Group 

Method 

Physics 

Type 

Grid 

Cont. Vol. 

1 

A 

CFL3D 

Euler 

Struct 

SI 

10,141,696 

2 

A 

CFL3D 

Euler 

Struct 

S2 

1,267,712 

3 

A 

CFL3D 

Euler 

Struct 

S3 

158,464 

4 

A 

CFL3D 

Euler 

Struct 

S4 

19,808 

5 

A 

FUN3D 

Euler 

Mixed 

M100 

5,363,065 

6 

A 

FUN3D 

Euler 

Mixed 

M125 

2,742,341 

7 

A 

FUN3D 

Euler 

Mixed 

M156 

1,411,515 

8 

A 

FUN3D 

Euler 

Mixed 

M200 

695,709 

9 

A 

FUN3D 

Euler 

Tet 

T100 

5,363,065 

10 

A 

FUN3D 

Euler 

Tet 

T125 

2,742,341 

11 

A 

FUN3D 

Euler 

Tet 

T156 

1,411,515 

12 

A 

FUN3D 

Euler 

Tet 

T200 

695,709 

13 

B 

FUN3D 

Euler 

Tet 

VGRID 

16,700,946 

14* 

C 

USM3D 

Euler 

Tet 

T100 

31,801,283 

15 

C 

USM3D 

Euler 

Tet 

T125 

16,191,339 

16 

C 

USM3D 

Euler 

Tet 

T156 

8,288,192 

17* 

C 

USM3D 

Euler 

Tet 

T200 

4,056,575 

18 

C 

USM3D 

Euler 

Tet 

PW/MCAP 

6,033,621 

19 

C 

USM3D 

SA 

Tet 

PW/MCAP 

7,864,672 

20* 

D 

WOLF 

Euler 

Tet 

Adapted 

759,048 

21* 

D 

WOLF 

Euler 

Tet 

Adapted 

860,010 

22* 

D 

WOLF 

Euler 

Tet 

T100 

5,363,065 

23* 

D 

WOLF 

Euler 

Tet 

T125 

2,742,341 

24* 

D 

WOLF 

Euler 

Tet 

T156 

1,411,515 

25* 

D 

WOLF 

Euler 

Tet 

T200 

695,709 

26 

E 

Eugenie 

Euler 

Tet 

T100 

5,363,065 

27 

E 

Eugenie 

Euler 

Tet 

T125 

2,742,341 

28 

E 

Eugenie 

Euler 

Tet 

T156 

1,411,515 

29 

E 

Eugenie 

Euler 

Tet 

T200 

695,709 

30 

F 

CEDRE 

Euler 

Tet 

T100 

31,801,283 

31 

F 

CEDRE 

Euler 

Tet 

T200 

4,056,575 

32 

F 

elsA 

Euler 

Struct 

SI 

10,141,696 

33 

G 

FAS IP 

Euler 

Struct 

Gridgen 

9,848,768 

34 

H 

SU2 

Euler 

Tet 

Output 

229,000 

35 

H 

SU2 

Euler 

Tet 

Adapted 

1,197,548 

36 

H 

SU2 

Euler 

Mixed 

M100 

5,363,065 

37 

H 

SU2 

Euler 

Mixed 

M125 

2,742,341 

38 

H 

SU2 

Euler 

Mixed 

M156 

1,411,515 

39 

H 

SU2 

Euler 

Mixed 

M200 

695,709 

40 

J 

UPACS 

Euler 

Struct 

SI 

10,141,696 

41 

J 

UPACS 

Euler 

Overset 

Fixed 

27,646,000 

42 

J 

TAS 

Euler 

Mixed 

M100 

22,423,283 

43 

J 

TAS 

Euler 

Mixed 

M125 

11,451,499 

44 

J 

TAS 

Euler 

Mixed 

M156 

5,877,312 

45 

J 

TAS 

Euler 

Mixed 

M200 

2,830,575 

46 

K 

ADCS 

Euler 

Struct 

Gridgen 

10,640,580 

47 

L 

FaSTAR 

Euler 

Mixed 

M100 

5,363,065 

48 

L 

FaSTAR 

Euler 

Mixed 

M125 

2,742,341 

49 

1 

FaSTAR 

Euler 

Mixed 

M156 

1,411,515 

50 

L 

FaSTAR 

Euler 

Mixed 

M200 

695,709 

51 

M 

Cart3D 

Euler 

Cart 

Output 

15,885,435 

52 

N 

Cart3D 

Euler 

Cart 

Adapted 

49,000,000 

53 

N 

Cart3D 

Euler 

Cart 

Fixed 

75,000,000 

54 

N 

CFD++ 

Euler 

Mixed 

M100 

22,423,283 

55 

O 

CartOver 

Euler 

Hybrid 

Fixed 

108,000,000 

56 

P 

Cart3D 

Euler 

Cart 

Fixed 

54,466,264 

57 

Q 

CFD++ 

Euler 

Mixed 

M100 

22,423,283 

58 

R 

LAVA 

SA 

Overset 

Fixed 

21,300,000 

59 

R 

LAVA 

SA 

Mixed 

Fixed 

12,100,000 

60* 

S 

OVERFLOW 

SA 

Overset 

Fixed 

29,600,000 


Updated or submitted after the workshop. 
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extraction location after traveling a longer distance than the centerline disturbances. This extra distance 
results in wave superposition at a more aft location that alters the signature shape. This signature change 
with extraction distance is known as near-field interference. When signatures are extracted sufficiently far 
from the model, a far-held propagation technique (sBOOM) is usable without correction because the acoustic 
sources appear to come from a single line. 

Extracted CFD pressure signatures were requested for a complete cylinder aligned with freestream how 
direction to enable a far-held correction methodology that accounts for a spanwise distribution of acoustic 
disturbances. 30 However, ground signatures and ground loudness in this summary apply sBOOM directly 
to data extracted at distances where the near-held interference of spanwise disturbances is signihcant. The 
far-held correction is omitted because it eliminates some of the differences between submissions. 

The scaled LM1021 wind tunnel model produces a louder PL than the hight vehicle it represents. The 
blade sting and the thickened boundary layers due to the lower wind tunnel Reynolds number distort the 
signature that is designed to be a quiet for the free-hight model at hight Reynolds number. CFD calculations 
on the full-scale free-hight model with far-held correction yield the PL in Table 6, which is quieter than the 
PL shown in this summary for the wind tunnel model without correction. The PL from two extraction 


distances varies by less than 

a tenth of a 

dB when 

near-field corrections are 

applied. 



Table 6. 

Full-scale LM1021 far-field corrected CFD. 



R/(b/2) extraction location 

PL at 4> = 0 

PL at <t> = 10 

PL at 0 = 20 PL at 0 = 30 

PL at 0 = 40 

PL at <£ = 50 

5.5 

85.04 

83.70 

84.62 85.67 

86.54 

86.26 

10.9 

85.00 

83.67 

84.52 85.52 

86.46 

86.28 


The 1/125 scale model has a reference length of L = 22.4 in. and a half span of 5/2 = 4.03 in. This 
yields a reference length to half span ratio of 5.6. The angle of attack is specihed as 2.1°. The provided 
grid contains geometry that was rotated to the specihed incidence, so it should be run at 0° angle of attack. 
The workshop committee requested turbulent calculations at a Reynolds number per inch of 212,000. This 
request is slightly lower than the actual test Reynolds number per inch of 375,000. This difference in Reynolds 
number had a negligible effect on the near-body signatures. 28 The increase in Reynolds number resulted in 
improved measurement accuracy and productivity. 28 

The participants provided submissions containing entire half cylinders, where pressure signatures are 
extracted every 10° at hve radial locations, Table 7. The submissions to this optional case are listed in 
Table 8. See section V.A for a description of the abbreviations that are common for the SEEB-ALR and 
LM1021 tables. 


Table 7. LM1021 pressure signature extraction locations. 


H 

H/L 

H/(b/ 2) 

<t> 



19.7 in. 

1.1 

4.9 

every 10° 

from 0° 

to 180° 

31.8 in. 

1.8 

7.9 

every 10° 

from 0° 

to 180° 

42.0 in. 

2.3 

10.4 

every 10° 

from 0° 

to 180° 

70.0 in. 

3.9 

17.4 

every 10° 

from 0° 

to 180° 

100.7 in. 

5.6 

25.0 

every 10° 

from 0° 

to 180° 


Submission 

Group 

Table 8. 

Method 

LM1021 submissions. 

Physics Type 

Grid 

Cont. Vol. 

1 

A 

FUN3D 

SA 

Tet 

WS 

2,395,158 

2 

B 

FUN3D 

SA 

Tet 

VGRID 

27,000,000 

3 

C 

USM3D 

Laminar 

Tet 

WS 

13,737,358 

4 

C 

USM3D 

SA 

Tet 

WS 

13,737,358 

5 

C 

USM3D 

SST 

Tet 

WS 

13,737,358 

6 

C 

USM3D 

SA 

Tet 

PW/MCAP 

65,400,000 

7* 

E 

Eugenie 

Euler 

Tet 

Fixed 

13,737,358 

8 

M 

Cart3D 

Euler 

Cut Cell 

Output 

57,000,000 

9 

O 

CartOver 

Euler 

Hybrid 

Fixed 

51,018,868 

10 

Q 

CFD++ 

SA 

Tet 

WS 

13,737,358 

11* 

Q 

CFD++ 

SA 

Tet 

Inflate 

23,463,802 


* Updated or submitted after the workshop. 
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VI. Near-Field Wind Tunnel Measurements 


Spatial and temporal variation in wind tunnel test section Mach, static pressure, and flow direction 
complicates near- field pressure measurements, because these noise sources are of the same magnitude as 
the signals the test is intended to measure. These disturbances have been known for decades 31 and are 
now better understood with the availability of time 32 and space 33 resolved measurements. Typical spatial 
disturbances are clearly illustrated by distortion in measurements of the flattop portion of the SEEB-ALR 
signature. 34 

A better understanding of the disturbances present in all supersonic wind tunnels 35 has led to a spatial 
averaging technique 34, 36 used in conjunction with a pressure measurement rail. 37 The evaluation and de- 
velopment of different rail configurations was also used to quantify the effect of humidity on wind tunnel 
measurements. 38 

The experimental signatures in this summary were the result of the spatial averaging technique, which 
provides mean signatures and variation among the ensemble of measurements used to compute the mean 
signature. This variation among the ensemble of measurements is expressed as a standard deviation at each 
x location. This spatial averaging technique was complicated by spatial variation in flow angularity 39 that 
caused variation in angle of attack and lift as the model was translated. 28 Bretl and Walker 40 propagated 
uncertainty derived from the near-field measurement variation and estimated uncertainties in the range of 
0.5 to 4 dB uncertainty in PL, depending on how the near-field variation was represented. 

Wind tunnel models for off-body near-field pressure measurements are generally smaller than force and 
moment models. These smaller models can amplify Reynolds number sensitivity and require characterization 
of the boundary layer state to increase the fidelity of simulation. 6 This is critical for models with small details 
(e.g. nacelle boundary layer diverters, 12 nacelles in the wake of wind tunnel mounting hardware 6 ). Boundary 
layer trip devices further complicate measurement and simulation because they are directly observed in near- 
body measurements and affect boundary layer state, which impacts aft-body signatures. 41,42 

A summary of the experimental measurements used in the workshop was provided by Cliff et al. 28 
Morgenstern, Buonanno, and Marconi 6 provided additional details of the LM1021 test. Makino and Noguchi 5 
performed measurements on an independently constructed DWB model. 

VII. SEEB-ALR Analysis 

The analysis for the SEEB begins by examining the submissions at the closest extraction distance and 
these signatures propagated to the ground. The variation in PL and validation metric measures is exam- 
ined with grid refinement and the statistics of these measures is examined. Finally, signatures that have 
statistically significant measures are examined. This same analysis process will be repeated for the DWB. 

VII. A. SEEB-ALR Near-Body Signatures 

The submitted signatures from Table 3 are shown as red lines in Fig. 8. The wind tunnel data has been 
shifted in x to align the bow shock to the middle of the participants. The mean of two wind tunnel test series 
are shown as the solid black lines. One standard deviation of the measurement ensemble variation is shown 
as the dashed black lines for each test series. All 64 submissions are shown in Fig. 8(a) and 30 submissions 
(excluding the coarser grids of the uniform refinement study) are shown in Fig. 8(b). The positions of the 
secondary nose shock, expansion, and the tail shock are predicted by the majority of the participants. The 
flattop pressure level is predicted within a standard deviation of the wind tunnel measurements, 28 < x < 37. 
The participant pressure signatures in this flattop region are at or above the lower pressure wind tunnel mean. 
The submissions are lower than the wind tunnel mean expansion pressures and slightly below one standard 
deviation of the wind tunnel measurements, 39 < x < 43. Eliminating the coarser uniformly refined grids 
does not appear to reduce the scatter of the results, Fig. 8(b). 

Some of the signatures provided do not have a point of inflection at the secondary nose shock at x = 28, 
Fig. 9. Participants C and D also provided signatures on the workshop provided tetrahedral grids that 
exhibited a secondary nose peak (not shown). Participant B only provided 15B, which used a custom grid 
and did not show the secondary shock. This missing secondary nose peak is likely due to geometry or 
grid differences, because the same CFD method using different grids showed the secondary peak. A possible 
explanation is the use of as-designed geometry which has as slightly different nose shape. 28 As-built geometry 
was provided for the workshop. 
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(b) Fine-grid submissions, excluding coarser uniform grids. 


Figure 8. 


SEEB-ALR near-body pressure signatures at H = 21.2 in. 
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Figure 9. SEEB-ALR, submissions exhibiting a single peak at H = 21.2 in. 


Submissions based on Euler calculations were requested for this case, but some participants submitted 
signatures extracted from turbulent RANS solutions. Figure 10 shows the turbulent submissions for the 
SEEB-ALR. The Spalart-Allmaras (SA) turbulence model was used for the 23C, 59Q, 61R, and 64S sub- 
missions. The Mentor Shear Stress Transport (SST) turbulence model was used for the 63R submission. 
These signatures (except for submission 23C) showed the secondary nose shock. As a group, the turbulent 
signatures had a slightly higher pressure in the aft expansion (closer to the mean wind tunnel measurements) 
than the Euler signatures, 39 < x < 43. 

VII. B. SEEB-ALR Ground Signatures 

All of the near-body SEEB submissions and the mean wind tunnel measurements at H=21.2 in. are propa- 
gated to the ground, Fig. 11(a). The propagation is performed with the method described in section II. B. 
The time axis of the ground signatures have been shifted so that all the initial pressure rises occur at 100 ms. 
The two mean wind tunnel measurement have slightly different flattop pressure levels, but nearly identical 
expansions and aft shock positions. The ensemble of submitted signatures produce ground signatures that 
are at or higher than the wind tunnel mean for the flattop and lower than the wind tunnel mean for the 
expansion. This is the same trend that is seen in the near-body signatures, Fig. 8(a). 

The coarser submissions in uniform grid refinement families are excluded in Fig. 11(b). The ensemble 
pressure scatter is not dramatically reduced by excluding the coarser grids. The is the same trend exhibited 
by near-body signatures in Fig. 8(b). 

VII. C. SEEB-ALR Uniform Refinement 

Ten sets of uniformly refined grids from eight of the nineteen participants were submitted for the SEEB-ALR. 
The participants were instructed to perform a calculation on the finest committee provided grid that they 
had the ability to run, and as many coarser grids in that family as possible. Some participants were not able 
to comply due to time constraints, and some solvers were not compatible with the provided grid types (e.g., 
Cartesian cut-cell methods, multiblock structured methods that can not handle singularities). 
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Figure 10. SEEB-ALR turbulent CFD submissions at H = 21.2 in. 


The PL of the propagated ground signatures is shown in Fig. 12(a). PL is calculated with the method 
described in section II. A. The characteristic grid spacing h is estimated as the number of control volumes 
to the negative one third power. The unstructured grids have an approximate doubling of control volumes 
between grids in the family. This results in successive unstructured grids with an approximate cube root 
of a half factor between successive h values. The structured grid has a factor of eight increase in control 
volumes (with a halving of h). Extrapolating the submission series to h = 0 (infinitely refined grid) yields 
approximately 91.7 PL (dB). There are no obvious trends for structured, mixed-element, or purely-tetrahedral 
grid submissions. 

The validation metric V is a measure of the difference between the submitted signatures and the mean 
wind tunnel measurement, Fig. 12(b). It is calculated with the method described in section II. C. This metric 
is expected to trend to a non-zero value with grid refinement because of the inherent differences between the 
mean wind tunnel measurement and submitted signatures (e.g. incomplete modeling of the physics with the 
Euler equations, rounding of the wind tunnel measurements due to averaging). Each SEEB-ALR uniform 
refinement series extrapolates to V ~ 0.0006 at h = 0 (infinitely refined grid). 

VII. D. SEEB-ALR Statistics 

The statistical approach of section III is used to compute the median and limits of the entire population of 
PL propagated from H = 21.1 in. in Fig. 13(a), which includes coarser uniformly refined grids by the same 
participant. The median is near 91.7 PL (dB), which agrees with the uniform grid refinement extrapolation. 
The two diamond symbols are the PL of both wind tunnel measurement means. The wind tunnel PL values 
are both slightly lower than the mean of the submissions, but within the limits of the participant statistics. 
All outliers have a larger PL than the upper limit. 

The fine grid submissions are shown Fig. 13(b) by removing the the coarser grids in the uniform grid 
refinement study from Fig. 13(a). The median and limit statistics are very similar for this smaller (but more 
independent) sample of the population. The outliers from participant E in Fig. 13(a) are on coarser grids; 
the finest grid submitted by participant E in Fig. 13(b) is not statistically significant among other fine-grid 
results. Signatures from the 21C and 431 submissions have a louder PL than the limits. Details of these near 
body signatures are examined in section VILE. 


16 of 43 


American Institute of Aeronautics and Astronautics 


Pressure (PSF) Pressure (PSF) 


0.8 r 



_n 8 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 i I 1 1 1 1 I 1 1 1 1 I 

‘ 50 100 150 200 250 300 350 

Time (ms) 


(a) All submissions. 



(b) Fine-grid submissions, excluding coarser uniform grids. 


Figure 11. SEEB-ALR propagated ground signatures from H = 21.2 in. 
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Figure 13. SEEB-ALR PL statistics propagated from H = 21.2 in. 
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The PL from both of the near body signature extraction locations are shown in Fig. 14(a). The wind 
tunnel PL is only propagated from H = 21.2 in. This slender, nominally axisymmetric body is expected to 
have very similar propagated PL values for these locations at 1.2 and 2.4 body lengths. Significant differences 
in PL from these two locations would indicate changes to the signature within the CFD domain. 

The differences between PL at H — 42.0 in. and H = 21.2 in. are shown in Fig. 14(b). All the outliers 
have a lower PL at the larger extraction location. Submission 25D is the only submission not involving a 
coarser grid in a uniform grid family. 
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(a) PL propagated from H = 21.2 in. and H = 42.0 in. 


(b) PL difference between H = 21.2 in. and H = 42.0 in. 
4A, 25D, 28D, 29D, and 52L are below the limits. 


Figure 14. SEEB-ALR PL statistics propagated from two near-body distances. 


VII. E. SEEB-ALR Outliers 

Submissions 21C and 431 near-field signatures are shown in Fig. 15. These submissions have a louder PL than 
the limits in Fig. 13(b). The 21 C Euler custom grid submission has a higher flattop pressure than the wind 
tunnel mean and signatures by the same participant using the same code on the committee-provide grids. 
Submission 431 is a linear method. This submission exhibits many of the small shocks in the measurement 
mean flattop region, but has a unique expansion shape, which exibits a much stronger expansion and recovery 
slightly before the aft shock. 

Submission 25D is the only submission (not a coarse grid in a uniform refinement family) that is outside 
of the limits in Fig. 14(b). The 25D near-body pressure signatures at H = 21.2 in. and H = 42.0 in. are 
scaled by yjH/L and shifted in x by the iL/3, Fig. 16. After this transformation, the two signatures should 
over plot because both extraction locations are sufficiently far from this axisymmetric body. However, the 
transformed H = 42.0 in. signature has a lower pressure value than the transformed H = 21.2 in. signature. 

VIII. DWB Analysis 

The same analysis steps applied to the SEEB-ALR are now applied to the second required workshop 
configuration, the DWB. It is slightly more complex than the SEEB-ALR and provides a slight variation in 
on- and off-track near-field signature because the wing has different impact on the signature at different </> 
locations. 

VIII. A. DWB Near-Body Signatures 

The submitted signatures from Table 5 are shown as red lines in Fig. 17 for H = 24.8 in. and <p = 0° 
(centerline). The mean of wind tunnel measurement ensemble is shown as the solid black line. One standard 
deviation of measurement ensemble variation is shown as the dashed black line. The wind tunnel data has 
been shifted in x to align the strong expansion to the median of the participants. All 60 submissions are 
shown in Fig. 17(a) and 31 submissions (excluding the coarser grids of the uniform refinement study) are 
shown in Fig. 17(b). Without the submission on coarser grids, the positions of the shocks and expansions 


19 of 43 


American Institute of Aeronautics and Astronautics 


(dp/pinf)*sqrt(H/L) 



X 


Figure 15. SEEB-ALR fine-grid PL outliers at H = 21.2 in. 


0.010 r 


0.005 - 


0.000 - 


-0.005 - 


- 0.010 - 


- 0.01 


- 0.2 



j I i i i I i i i I i i i I i i i I i i i I i i i I 

0 0.2 0.4 0.6 0.8 1 1.2 

(X-H*Beta)/L 


Figure 16. SEEB-ALR fine-grid PL difference outlier. 
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are predicted by the majority of the participants. The first two shocks are due to the nose and the wing. 
The series of small shocks starting at X = 1.08 m are due to the expansion and re-compression of the base 
of the model. The pressure wave centered at X = 1.18 m, only seen on coarse grids of Fig. 17(a), is due to 
poor wind tunnel sting discretization on the coarser committee provided grids. 27 There is a large spread on 
the predicted level of overshoots and undershoots for the shocks. There is a small expansion in the middle 
of the aft wing shock that is predicted by the participants at X = 1.05 m but not shown in the mean wind 
tunnel measurements. 

VIII. B. DWB Ground Signatures 

The submitted signatures are propagated to the ground from H = 24.8 in. and 0 = 0° in Fig. 17. All 
submissions are shown in Fig. 18(a) and the coarser grids of the uniform refinement study are excluded in 
Fig. 18(b). The propagated submissions are shown in red and the propagated wind tunnel median is shown 
in black. The time axis of the ground signatures have been shifted so that all the initial pressure rises occur 
at 100 ms. The large variation in overshoots and undershoots and the complex signature due the the base 
flow seen in the near field have been filtered by the propagation process. The location of the middle and aft 
shocks shows more variation than the expansion slopes. This change in shock positions results in different 
minimum and maximum pressures of the main expansion. 

VIII. C. DWB Uniform Refinement 

Ten sets of uniformly refined grid submissions from eight participants are available for the DWB. The PL of 
the propagated ground signatures is shown in Fig. 19. The characteristic grid spacing h is estimated as the 
number of control volumes to the negative one third power. The DWB, Fig. 19, shows more variation than 
the SEEB-ALR, Fig. 12(a), when extrapolating the finest two grids from each series to zero h (infinitely 
refined grid). This extrapolation to zero h results in PL values between 95.2 and 95.6 dB. 

The grid convergence of the validation metric is examined in Fig. 20. The metric for the entire signature 
at H = 24.8 in. and 0 = 0° is shown in Fig. 20(a). When each V series is extrapolated to zero h (infinitely 
refined grid) there is less agreement of the extrapolated V than than SEEB-ALR V or DWB PL. To isolate 
the source of the larger spread in extrapolated values, the validation metric integral, Eq. (1), extent is limited 
to the forward portion of the signature to isolate it from the the impact of the base flow region. The whole 
extent is X = [0.8, 1.25] and the forward extent is X = [0.8, 1.0]. Restricting V to the forward portion of the 
signature did not reduce the scatter of the results, indicating that the base flow is not the only contributor to 
the scatter. The overshoots and undershoots of the strong shocks may be adversely impacting convergence 
of V with grid refinement. 

VIII. D. DWB Statistics 

The DWB PL propagated from H = 24.8 in. is shown for the centerline (Fig. 21), 30° off-track (Fig. 22), 60° 
off-track (Fig. 23), and 90° off-track (Fig. 24). The DWB sBOOM atmospheric propagation is performed 
directly down for each extracted signature, so every signature experiences the same atmospheric model. This 
approach is not the process commonly used to compute boom “carpets” where each signature is propagated 
at the same oblique angle where it is extracted. Using the same atmospheric model for all DWB signatures 
facilitates comparison between near-field signatures extracted at different 0-angles, including 90° off-track 
which would not reach the ground in a standard atmosphere. The sample median (solid black lines) and the 
propagated mean wind tunnel measurement (black diamonds) compare favorably for all 0. Both the median 
participant and experimental PL values and increase with 0. This is because the wing shock strength and 
spatial separation between the trailing edge wing shock and the base region shocks change with 0. 

All submissions are plotted in each (a) sub-figure and the coarse grids of the uniform refinement series 
are excluded from each (b) sub-figure. The sample median (solid lines) is insensitive to excluding the coarser 
refined grids, but the limits (dashed lines) grow slightly by approximately a tenth of a dB. This increase 
in the limits does not eliminate any fine-grid submissions that are outliers of all submissions. All DWB 
outliers are below the limits. Submission 4 A is the coarsest DWB grid in the workshop with 19.8 thousand 
cells/control volumes. Submission 33G is an outlier for all four 0-angles. Participant G indicates in their 
workshop presentation that the DWB structured grid wing tip singularity is an issue that negatively impacted 
the CFD tool. Submission 53N is an outlier at and above 0 = 30°. This submission uses a manually refined, 
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(b) Fine-grid submissions, excluding coarser uniform grids. 


Figure IT. DWB near-body pressure signatures at H = 24.8 and cf> = 0 °. 


22 of 43 


American Institute of Aeronautics and Astronautics 


Pressure (PSF) Pressure (PSF) 



(a) All submissions. 



(b) Fine-grid submissions, excluding coarser uniform grids. 

Figure 18. DWB propagated ground signatures from H = 24.8 and 4 > — 0 °. 
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Figure 19. DWB PL uniform grid refinement propagated from H = 24.8 and 0 = 0 °. 
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(a) The entire signature. (b) The forward portion of the signature. 

Figure 20. DWB validation metric uniform grid refinement at II = 24.8 and 0 = 0 °. 
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(b) Fine-grid submissions; 34H is below the limits. 


Figure 21. DWB PL statistics propagated from H = 24.8 in. and 0 = 0°. 
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(a) All submissions; 4A, 34H, and 53N are below the limits. 


(b) Fine-grid submissions; 34H and 53N are below the lim- 
its. 


Figure 22. DWB PL statistics propagated from H = 24.8 in. and </> = 30°. 
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(a) All submissions; 4A, 33G, 34H, and 53N are below the 
limits. 


(b) Fine-grid submissions; 33G, 34H, and 53N are below 
the limits. 


Figure 23. DWB PL statistics propagated from H = 24.8 in. and cj) = 60°. 
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(a) All submissions; 33G, 34H, and 53N are below the lim- 
its. 


(b) Fine-grid submissions; 33G, 34H, and 53N are below 
the limits. 


Figure 24. DWB PL statistics propagated from H = 24.8 in. and cj) = 90°. 
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centerline shock aligned Cartesian grid that has low resolution of the DWB base flow. The output adapted 
34H submission is an outlier at and above 0 = 60° and is the coarsest grid (229,000 control volumes) of the 
fine-grid submissions. Other participant H submissions based on a larger adapted and workshop provided 
mixed-element grids are not statistically significant at any 0. 

The PL of propagated submissions from three centerline locations is shown in Fig. 25(a) with PL of 
propagated mean wind tunnel measurements. The winged DWB configuration is insensitive to extraction 
distance over this range of extraction distances (3.1, 3.6, and 4.6 body lengths or 15.6, 18.3 and 23.4 half 
spans) at this non-lifting flow condition. The difference between the PL propagated from H = 31.8 in. and 
H = 24.8 in. is shown in Fig. 25(b). The scatter of the DWB submissions is lower than the scatter of the 
SEEB-ALR submissions, Fig. 14(b). The ratio of extraction distances is also slightly less. The median of the 
PL differences is similar to the mean wind tunnel measurement PL difference. The coarse-grid 4A submission 
is above the limits and the 51M submission is slightly below the limits. 
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(a) PL propagated from three near-body distances. 


(b) Difference between H = 31.8 in. and H = 24.8 in. 4A 
is above and 51M is below the limits. 


Figure 25. DWB PL statistics propagated from three near-body distances on the centerline. 


VIII. E. DWB Outliers 

Submissions 33G, 34H, and 53N are identified as statistical fine- grid outliers in PL at multiple 0. These 
submissions are not all outliers at every 0-angle, but are shown at all 0 to illustrate 0-trends. The near- 
field CFD pressure signatures are shown on the centerline (Fig. 26), 30° off-track (Fig. 27), 60° off-track 
(Fig. 28), and 90° off-track (Fig. 29) with mean wind tunnel measurements. The wind tunnel mean shows 
the separation of the wing and base flow shocks and the steepening of the wing expansion with increasing 
off-track angle. Submission 33G has a unique shock in the middle of the main expansion that grows with 
0. Participant G indicated that the wing tip structured grid singularity may have adversely affected the 
CFD solver. The adapted-grid 34H submission predicts a lower pressure in the expansions than the wind 
tunnel mean, which has a constant expansion slope. Submission 53N lacks distinguishing features that would 
indicate why it is an outlier. 

One outlier for the difference in PL at two extraction distances remains after coarser uniformly refined 
grids are excluded. The near- field signature for the remaining outlier 51M is shown in Fig. 30. The slope 
of the main expansion is less steep for the 51M near-field centerline signature at H = 31.8 in. Also, the aft 
wing and base flow shocks have merged. This is the same behavior exhibited by most participants. It is 
unclear what subtle difference results in the reduction of 0.2 PL (dB) from H = 24.8 in. to H — 31.8 in. 

IX. LM1021 Analysis 

Some of the analysis steps from the SEEB-ALR and DWB are applied to the LM1021. The LM1021 has 
fewer submissions, so the previously used statistical methods are not be applied to the LM1021. The LM1021 
submissions show a larger amount of variation in near-field and ground signatures than the SEEB-ALR or 
DWB, which may be due to the complexity of the model or the tailored nature of the low-boom signature. 
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Mgure 26. DWB fine-grid PL outliers: 


33G, 34H, and 53N at II = 24.8 in. and 0 = 00 
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Figure 27. DWB fine-grid PL outliers: 33G, 34H, and 53N at H = 24.8 in. and 0 = 30 
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Mgure 28. DWB fine-grid PL outliers: 


33G, 34H, and 53N at II = 24.8 in. and 0 = 60 
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Figure 29. DWB fine-grid PL outliers: 33G, 34H, and 53N at H = 24.8 in. and 0 = 90 
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Figure 30. DWB fine-grid PL difference outlier: three near-body distances for scaled 51M. 


The variation near-field is also reflected in a large variation in PL and CSEL. 

IX. A. LM1021 Near-Body Signatures 

To introduce the LM1021, all of the submitted near-field signatures on the centerline at H = 31.8 in. are 
shown in Fig. 31 as gray lines; CFD was requested at 2.10° angle of attack. Mean wind tunnel measurements 
are shown from a number of test series in Fig. 31 for conditions near H = 31.8, each with a different mean 
angle of attack. The 2.30° angle of attack mean signature is computed with amplitude scaling, age scaling, 
reverse aging, auto reference clipping, and other processing features. Cliff et al. 28 describes the method used 
to compute the other mean signatures. The wind tunnel mean pressure level at the forward portion of the 
fuselage ( X = [3.3, 3.9]), the forward portion of the wing ( X = [3.9, 4.4]), the wing shock position ( X = 4.5), 
the nacelles ( X = [4.7, 4.9]), and the tail ( X > 5) vary with angle of attack, wind tunnel test series, and 
spatial averaging technique. 34,36 The angle of attack of the wind tunnel test is corrected with an upright and 
inverted measurement of lift. 21 This process is complicated by flow angularity variation in the wind tunnel 
test section that results in 0.16° peak-to-peak variation in angle of attack 28 as the model is translated to 
perform spatial averaging. The 2.10° angle of attack specified for workshop CFD computations is based on 
CFD performed before the workshop. 6 These 2.10° angle of attack pre- workshop simulations compared well 
to 2.30° angle of attack wind tunnel measurements. There are also data reduction updates that indicated 
that the wind tunnel angle of attack may be lower than initially estimated. 

All the submitted signatures from Table 8 are shown in Fig. 32 at H = 31.8 in. on the centerline. There 
is a significant variation in pressure level over the entire signature and in the details of the main shock and 
expansion in the influence region of the nacelle ( X = [4.6 — 4.9]). The variation in pressure level is largest at 
the tail of the configuration, especially the pressure recovery behind the model ( X > 5.3). The differences in 
pressure level ahead of the wing ( X < 4) may be due to differences in boundary layer displacement thickness. 
The submissions are grouped by physical model in subsequent figures. 

The laminar and turbulent submissions are shown in Fig. 33. The Euler submissions are shown in Fig. 34. 
Referring to all submissions in Fig. 32, the laminar and turbulent submissions have a higher pressure for 
the forward portion of the signature than the Euler submissions. The predicted laminar and turbulent main 
shock location is ahead of the Euler main shock location. 
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Figure 32. All LM1021 near-field submissions at H = 31.8 in., cf> = 0°. 
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There is little variation between the three Euler signatures (Fig. 34), except at the main shock, nacelle, 
and tail. There is a shock in the main expansion that is not seen in the laminar and turbulent submissions 
or wind tunnel mean measurements. This difference, unique to Euler calculations, may be due to the 
boundary-layer choking the wing- nacelle gap flow and preventing this shock from forming. 12 
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Figure 33. Laminar and turbulent LM1021 near-field submissions at H = 31.8 in., 0 = 0 °. 


Participant C provided three calculations on the same grid, Fig. 35. These include a laminar calculation 
(3C) and turbulent calculations with the SA (4C) and SST (5C) turbulence models. The two turbulent 
simulations produce equivalent signatures. The laminar simulation has slightly lower pressure for the forward 
portion of the signature X < 4.4 and the main expansion at X = 4.8. 


IX. B. LM1021 Ground Signatures 

All the near-field signatures from Fig. 32 are propagated to the ground in Fig. 36. The signatures are 
propagated as provided by the participants without the application of a far-held correction 30 that would 
be appropriate for signatures extracted this close to a vehicle with a spanwise distribution of acoustic 
disturbances. The signatures have been shifted in time to align the 0.2 PSF pressure rise to 100 ms. The 
spread in forebody pressure seen in Fig. 32 results in a large spread in peak pressure in the first 50 ms of 
the ground signature. There is a large variation in the center shock location on the ground. The variation is 
smaller in the expansion after this shock. The location and strength of the final shock shows a large amount 
of variation. 

The laminar and turbulent submissions shown in Fig. 33 are propagated to the ground in Fig. 37. The 
three Euler submissions shown in Fig. 34 are propagated to the ground in Fig. 38. The peak pressure of the 
laminar and turbulent submissions is higher and shows more variation than the three Euler submissions. The 
differences in the near-field main shock and nacelle region have manifested themselves as differences in the 
middle ground shock locations, where the middle Euler ground shock locations are all aft of the laminar and 
turbulent middle ground shock locations. With the exception of submission 2B, the laminar and turbulent 
tail shock locations are aft of the Euler tail shock locations. Submission 7E differs from the other two Euler 
signatures in shock locations and the shape of the final shock. 

The participant C laminar near-field signature (3C) and turbulent near-field signatures with the SA (4C) 
and SST (5C) turbulence models are shown in Fig. 35 for the near-field and propagated to the ground in 
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Figure 34. Euler LM1021 near-field submissions at II = 31.8 in., 0 = 0 °. 


3C 

4C 



X 


Figure 35. Laminar and turbulence model effects on LM1021 near-field signatures at H = 31.8 in., (f) 
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Figure 36. All LM1021 ground signatures from H = 31.8 in., 0 = 0 °. 
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Figure 37. Laminar and turbulent LM1021 ground signatures from H = 31.8 in., 0 = 0 °. 
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Figure 38. Euler LM1021 ground signatures from H = 31.8 in., 0 = 0°. 


Fig. 39. The two turbulence models only show differences in the nacelle influence region (t = 270 ms). The 
laminar signature differs from the turbulent signatures primarily in shock locations. The laminar submission 
3C shock located at t = 270 ms is in the most aft position of any submission. 

IX. C. LM1021 Loudness 

The PL and CSEL are computed from the propagated LM1021 ground signatures. In general terms, PL 
emphasizes the higher frequencies between 1,000 and 10,000 KHz and CSEL covers a lower and wider 
frequency range between 100 Hz and a 5,000 KHz. These measures are shown as a function of near-field 
extraction distance on the centerline in Fig. 40. These measures are also presented in Fig. 41 as a function 
of off-track angle (j) for the signatures extracted at R/(b/ 2) = 7.9. In contrast to the DWB, each off-track 
angle for the LM1021 is propagated obliquely through the atmosphere creating a “boom carpet.” 

The signatures are propagated as provided by the participants without the application of a far-field 
correction 30 that would be appropriate for signatures extracted this close to to the vehicle. This correction 
is omitted because it eliminates some of the differences between submissions that this summary intends to 
quantify. The combination of no far-field corrections, inclusion of blade model support, and simulating the 
wind tunnel Reynolds number results in a PL that is louder than a full-scale free-flight model, see section V.C 
for details. 

A series of uniformly refined grids were not provided for the LM1021 configuration. The lack of signatures 
computed on a series of uniformly refined grids combined with the small number of submissions hinders the 
extraction of patterns from the LM1021 loudness, but the following observations are offered. The viscous 
methods approach 91-94 dB PL and 99-97 dB CSEL at far R with the exception of 2B, 10Q, and 11Q. 
These submissions have a lower pressure than other viscous submissions in the front portion of the near- 
field signature (Fig. 33) that results in a lower peak pressure (Fig. 37). Submission 2B has a high pressure 
level in the aft part of the near-field signature that results in a forward tail shock location (near the tail 
shock location of the Euler methods). The Euler methods are near 86 dB PL and 95 dB CSEL with the 
exception of 7E. Submission 7E has a unique ground tail shock position and shape, Fig. 36 and Fig. 38. 
Submission 1A has the greatest increase in PL with R and submission 10Q has a unique drop in CSEL 
level with increasing R. With the exception of 10Q, all the participants show less R variation in CSEL than 
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Laminar and turbulence model effects on LM1021 ground signatures propagated from H = 31.8 in., 0 = 0 °. 
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(a) PL. (b) CSEL. 

Figure 40. LM1021 loudness variation with near-body distance at <f> = 0°. 
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(a) PL. (b) CSEL. 

Figure 41. LM1021 loudness variation with off-track angle at H = 31.8 in. 


PL. CSEL emphasizes more of the low frequency content of the ground signature than PL, which indicates 
the propagated ground signatures have more high frequency changes with extraction distance R than low 
frequency changes. Sensitivity analysis of LM1021 PL and CSEL measures indicates that both of these 
measures are sensitive to angle of attack (lift). 43 

The trends with off-track angle at R/(b/2) = 7.9 are shown in Fig. 41. The PL and CSEL of all 
submissions decrease between 0° and 20° off track, where the variation in PL is the least. The trend with 
increasing off-track angle is less consistent above 20° off track with one group growing louder and the other 
staying nearly constant, The same louder and quieter groups are present in CSEL. The PL differences between 
10Q and 11Q, seen on the centerline, trend toward zero at 50° off track. Submission 6C with a custom grid 
is quieter than other participant C submissions on the committee provided grid. A possible explanation of 
6C off-track quieting could be differences in the off-track grid resolution between the committee provided 
and 6C grid. 

The two submissions from participant Q are the same CFD method applied to the committee provided grid 
(10Q) and a custom grid (11Q). These two scaled and windowed near-field signatures are nearly identical 
except for aft pressure peaks, Fig. 42. These small differences result in very fine changes to the ground 
signatures, Fig. 43. However, there is approximately a 3 dB PL offset over a range of R/(b/ 2) extraction 
locations, Fig. 40. Examining these ground signatures in the frequency domain (Fig. 44) shows that they 
are identical below 100 Hz, but 11Q has more energy above 100 Hz. Applying an A- weighted digital filter to 
the ground signatures isolated the differences to the sharpness of the tail shock. 44 The A-weighted loudness 
shows the same 3 dB difference between the Q submissions as PL, but CSEL only shows a 0.1 dB difference 
between 10Q and 11Q at R/(b/ 2) = 7.9. This is due to CSEL including more lower frequency energy than 
PL. These lower frequencies are identical between the two ground signatures, Fig. 44. The subtle differences 
between 10Q and 11Q may be useful as an evaluation tool for propagation and loudness measures. 

IX. D. LM1021 Far-Field Correction 

The spanwise distribution of disturbance sources cause near-field interference because these disturbances 
reach a near- field extraction location after traveling different distances. The CFD extraction location should 
be many times the semi-span of the model to minimize this near-field interference effect. When signatures 
are extracted sufficiently far from the model, a far-held propagation technique (sBOOM) is usable without 
correction because the acoustic sources appear to come from a single line. 

CFD data is requested for an complete cylinder around the configuration to enable a far-held correction 
methodology that accounts for a spanwise distribution of acoustic disturbances. 30 The average difference 
in PL of signatures propagated with and without this correction is shown in Fig. 45(a) as a function of 
R/(b/ 2). The impulse of this near-held signature correction is shown in Fig. 45(b). Both measures of the 
far-held correction decay with R. Far-held correction impulse shows a very steady hyperbolic convergence. 
The convergence of PL is less regular because the far-held correction changes the shape of the signature 
at different R/(b/ 2) locations. This shape change may result in an increase, decrease, or no change to PL. 
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42. LM1021 full-scale near-field signatures for participant Q at H = 31.8 in., cf> = 0°. 



Figure 43. LM1021 ground signatures for participant Q from H = 31.8 in., (f) = 0°. 


37 of 43 


American Institute of Aeronautics and Astronautics 


120 


N 

X 

m 

"D 


100 - 


80 - 


c n 

o 60 

Q 

ro 


o 

CD 

a. 

CO 

>% 

o> 

k— 

CD 

c 

LU 



40 - 


20 - 


0 - 


-20 


10 u 


10Q 
1 1Q 


_i i i i 1 1 1 1 


_i i i i 1 1 1 1 


10 1 



Frequency, Hz 


Figure 44. LM1021 ground signature spectra for participant Q from H = 31.8 in., <f> = 0°. 


Intermediate extraction locations may be useful for quantifying design changes, but a large R/(2/b) > 25 
extraction distance is required for fax-field propagation methods to obtain loudness measures within a single 
dB without a far-field correction. 
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Figure 45. LM1021 far-field correction decay with near-body distance. 


IX. E. LM1021 Equivalent Area 

Integrating the CFD signature into an equivalent area by the method of Whitham 45 as extended by Walk- 
den 46 quantifies the overall integrated strength of each signature. Equivalent area is computed for each 
participant that provided an entire half cylinder of extracted pressures, Fig. 46. The equivalent area in- 
creases linearly with total lift (boundary layer displacement is a small additional contribution), so the CFD 
differences correspond to a ±15% variation in lift for these submissions. Submissions 7E and 90 are Euler, 
submission 3C is laminar, and the other submissions are viscous. These two Euler submissions have a lower 
equivalent area than the laminar or turbulent submissions. The average equivalent area of the viscous solu- 
tions is approximately 80 square feet more than the average Euler submission. The ±6% variation between 
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the viscous submissions may be due to differences in lift or boundary layer state. 



Figure 46. LM1021 equivalent area. 


The equivalent area calculation includes scaling for the square root of R dependency of pressure. Each 
submission should have constant area value when calculated from any of the 5 provided R/(b/2) distances. 
These solutions show mostly constant areas except for an abrupt reduction between 17.4 and 25 R/(b/ 2) 
for all but two of the solutions (11Q and 90). Submissions that exhibit a change in equivalent area may 
indicate that grid resolution is insufficient for accurate results at theses larger distances. The 7E and 90 
Euler submissions have the same equivalent areas at R/(b/2) <= 8, but very different PL and CSEL in 
Fig. 40. 


X. Recommendations 

The iterative convergence of supersonic flows is complicated by the use of reconstruction limiters. 47 
However difficult, an iterative convergence criteria should be established or iterative convergence histories 
should be requested from participants. Insufficient iterative convergence was detected in the workshop 
uniform grid refinement studies of some participants, which encouraged the participants further converge the 
CFD for this summary. The convergence of forces is often monitored as a surrogate, but participant J showed 
in their workshop presentation that near-field signatures are still converging after the forces have reached 
steady state. This indicates that monitoring the iterative convergence of the solution at the extraction 
location may be a useful and possibly necessary tool. 

Length units should be standardized for future workshops. The SEEB-ALR and LM1021 geometry and 
grids were provided in inches, but the DWB grid and geometry was provided in meters. The SEEB-ALR 
signatures were required in inches, the DWB signatures were required in meters, and the LM1021 signatures 
were required in feet. The use of different systems of units unduly confused participants and complicated 
the analysis of submissions. 

Having a common grid system that the most participants could run was critical to understanding the 
results provided on custom and adapted grids. These committee provided grids identified problems with 
iterative convergence, post processing, and grid adaptation that would have been impossible to isolate if the 
participants had only used custom grids. Approximately half of the participants ran multiple grids in the 
workshop required uniform grid refinement series. These uniformly refined signatures were invaluable for 
understanding the convergence of PL and V measures. A continued effort should be made to ensure that as 
many participants as possible can run multiple grids in a uniformly refined grid series with as many CFD 
methods as possible. LM1021 submissions on multiple uniformly refined grids may have aided the discovery of 
the issues that caused the large variation in loudness measures of signatures provided by different particpants 
on the same grid. 

A grid size recommendation for adapted grids should be established. Many adapted grid submissions 
utilized very small grids. This contributed to faster execution time, which is desirable in a production analysis 
or design setting. However, under resolved signatures due to a small number of control volumes could be 
misinterpreted as a failure of the adaptive grid scheme. Requesting both intermediate and final adapted 
signatures would indicate how the PL and V measures converged as compared to uniform grid refinement 
and would verify the correctness of the adaptive method. 12,48 > 4 9 

The LM1021 signature and loudness metrics are sensitive to lift. 43 The computed lift should have been 
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gathered from participants and reported as a tool to understand the differences in near-body signatures 
and propagated ground signature loudness. An alternative would have been to run the CFD simulations 
at a specified lift condition and request angle of attack, which was done in the DPW series. The LM1021 
extraction macro has about a quarter of resolution in terms of points per body length than the SEEB- 
ALR or DWB. The sensitivity of propagated signature to the resolution of the extraction macro should be 
understood. 

The far-held correction 30 was significant for the LM1021 at the extraction distances where high quality 
wind tunnel data was available. This indicates that directly propagating the participant submissions with an 
augmented Burgers code from these locations is not as appropriate as it is with the SEEB-ALR and DWB. 
The wind tunnel mounting hardware and the thicker boundary layer at the wind tunnel Reynolds number 
negatively impacted the signature of the model that is designed to be a quiet full-scale free-hight model. 
The differences between wind tunnel test and fight would be a relevant topic for future study and inclusion 
in a subsequent workshop. 


XI. Next Steps 

The NASA Supersonics Project sonic boom goal for future supersonic aircraft is 65 to 70 PL (dB). 50 
Henne 15 proposes a small quiet supersonic jet (QSJ) with 68 to 85 PL (dB). Kubota 51 suggests a target of 88 
PL (dB) in order to attain an acceptability of 90% in laboratory experiments. The configurations for SBPW- 
1 were significantly louder: the DWB was approximately 96 PL (dB), the SEEB-ALR was approximately 
92 PL (dB), and the LM1021 was 84-94 PL (dB). Examining submissions from quieter configurations with 
the methods developed in conjunction with SBPW-1 would maintain the relevancy of future sonic boom 
workshops. 

Having a range of geometric complexity was important to the participants of SBPW-1. A future workshop 
should include both an axisymmetric configuration and a full-configuration (wing, body, tail, and nacelle) 
with complex near-field signatures that produce mid-80s to mid-70s PL on the ground. Having axisymmetric 
configuration that produces a signature with the complexity of the LM1021 may help to isolate geometry 
complexity as a source of PL variation. A formation of multiple axisymmetric bodies may isolate the spanwise 
distribution of acoustic disturbances as a source of PL variation. 

Unfortunately, a cursory survey of publicly available geometries with wind tunnel measurements did not 
yield a sufficiently quiet configuration. This leaves the possibility of performing the next workshop on a new 
configuration. There is a risk that this new configuration may never be tested in a wind tunnel, but DPW-4 52 
accepted this risk by holding workshop with an untested model and then successfully advocated for a test 
of that model after the workshop. Model supports and thicker boundary layers at lower Reynolds number 
required for supersonic wind tunnel testing can impact a design that is optimized for full-scale free- flight. It 
may be desirable to study a vehicle in both full-scale free- flight and wind tunnel test configurations. 

The propagation of the near-field signatures to the ground and the loudness calculation from ground 
signatures can introduce additional uncertainty sources to sonic boom prediction. West et al. 43 described 
a method of propagating uncertainties in the near-field CFD and atmospheric boom propagation, which 
showed that uncertainty in the atmospheric conditions can contribute more to PL uncertainty than the near- 
field CFD alone. Therefore, future workshops should include an examination of propagation and loudness 
calculation methods to gain additional insight on the entire uncertainty budget of sonic boom prediction. 

The impact of undershoots and overshoots, commonly present in near-field CFD signatures, has not 
been completely quantified for sonic boom prediction. Propagation methods, particularly methods that 
include molecular relaxation, filter the high frequency content present in near-field signature undershoots 
and overshoots. Direct frequency analysis of the near-field signature could be used to quantify impact 
of high frequencies the how they are filtered by the propagation method. 53 Digital filtering of near-field 
and propagated signatures could be used understand what sections of the signature contribute most to the 
variation in high frequencies that could be impacting PL. 44 

This summary focused on the PL noise measure due to its published correlation to annoyance 17 and its 
specification as a design target. 15,50,51 Recent experiments have indicated that other measures are correlated 
with annoyance, particularly for sonic booms experienced indoors. 54-56 Future workshops should consider 
these measures in addition to PL to understand if they place different requirements on near-field CFD or 
propagation methods because these measures emphasize different portions of the frequency spectrum. 
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XII. Conclusions 


A summary and statistical analysis of the First AIAA Sonic Boom Prediction Workshop (SBPW-1) has 
been presented. Submissions from 19 groups are presented for two required and one optional configuration. 
These configurations captured a range of geometric complexity. They are an axisymmetric body with flat- 
top ground signature (SEEB-ALR), a simple non-lifting wing body with a N-wave signature (DWB), and a 
complex configuration with wing, tail, flow-through nacelles, and blade mount sting (LM1021). 

The participants provided near-field pressure signatures extracted from CFD computations at specific 
locations. The SEEB-ALR had 64 submission, the DWB had 60 submissions, and the optional LM1021 
had 11 submissions. Many participants generated custom grids and solution adapted grids based on their 
best practices. The majority of participant groups provided multiple submissions computed with a series of 
uniformly refined committee provided grids. Submissions produced on these uniformly refined grids enabled 
a unique opportunity to study the grid convergence of a wide range of CFD methods for metrics of interest 
to the sonic boom research community. 

The provided near-field signatures were compared to each other and wind tunnel measurements. A valida- 
tion metric that quantifies the difference between the provided signatures and wind tunnel measurements was 
computed. The near-field CFD and wind tunnel signatures were propagated to the ground for comparison 
and the loudness measure PL was computed. 

Statistics of validation metrics and loudness measures were computed for the two required SEEB-ALR and 
DWB configurations. The SEEB-ALR had a quieter median 91.8 PL (dB) than the DWB on the centerline, 
95.5 PL (dB). The DWB median increased with off-track angle. All the median values for the required 
cases were within 0.5 dB of propagated wind tunnel measurements and often within 0.1 dB. The inclusion 
or exclusion of the coarse grids of the uniform refinement study had a negligible impact on median and limit 
statistics. The limits were about ±0.8 dB for the SEEB-ALR and ±0.2 dB for the DWB. The uniform grid 
refinement study on the committee provided grids showed better convergence for the SEEB-ALR than the 
DWB for PL and the validation metric. Scatter was reduced with grid refinement for the configurations and 
metrics examined, except for the DWB validation metric. The DWB model produces strong shocks, which 
may have impeded the convergence of the validation metric. 

The optional LM1021 had eleven submissions with larger scatter in PL than the two required cases. The 
differences between Euler, laminar, and turbulent submissions were noted in the near-field and propagated 
ground signatures. The PL measure had a large variation between methods and at different extraction 
location for the same method. The CSEL measure had less variation with extraction location. The louder 
PL submissions on the centerline were also loudest off track. A pair of signatures was identified that appeared 
very similar at the near field, but had almost a 3 dB PL difference when propagated to the ground with 
sBOOM. 

The decay of the far-field correction with extraction distance was documented in terms of impulse and its 
impact on PL for this configuration. Impulse smoothly decreased with distance, but the PL correction had 
different behavior at each off track angle. The PL correction has more than one dB at some off-track angles at 
25 half spans. Equivalent area computed from near-field signatures was presented as a function of extraction 
distance. The Euler submissions had a lower equivalent area than viscous submissions. Most submissions 
exhibited a drop in area between the furthest two extraction locations indicating that grid resolution was 
insufficient to maintain the signature. 

Recommendations for future sonic boom workshops were provided with a discussion of the next steps in 
preparation for a potential follow-on workshop. These recommendations identified weaknesses in the analysis 
methods and opportunities to facilitate better understanding at future workshops. 
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